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ABSTRACT 


A study  of  non-isoenergetic  turbulent  jet  mixing  between 
two  streams  has  been  conducted.  Using  a previously  derived 
theoretical  analysis  for  ducted  mixing,  an  experimental  investi- 
gation was  performed  to  verify  this  theory  and  to  determine 
the  non-isoenergetic  turbulent  jet  mixing  characteristics  in 
a constant  area  duct.  Temperature  profiles  were  measured  at 
several  axial  locations  in  the  duct  for  both  a concentric  and 
an  eccentric  configuration.  It  was  determined  that  the 
theoretical  and  experimental  temperature  profiles  agreed  fairly 
well  for  both  cases,  although  the  concentric  case  showed  better 
agreement  than  the  eccentric  case.  It  was  also  determined  that 
a new  constant  of  turbulence  in  the  initial  region  was  needed 
for  non-isoenergetic  mixing,  mixing  is  generally  more  rapid  than 
the  theory  predicted,  the  initial  temperature  difference  between 
the  two  streams  did  not  have  much  effect  on  the  rate  of  mixing 
and  a higher  area  ratio  produced  better  agreement  between  the 
theory  and  the  experimental  data.  It  was  concluded  that  the 
theory  was  good  for  a fairly  simplified  analysis. 
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NOMENCLATURE 


Symbol 


A 

AR 

C 


e 

e 

L 

1 

M 


m 

m 

n 

P 

R 


r 

r 


AT 

U 


2 

Cross  sectional  area  of  the  duct,  ft 
Area  ratio 

Experimental  constant  of  turbulence  (main  region) 

Experimental  constant  of  turbulence  (initial  region) 

Specific  heat  at  constant  pressure,  Btu/lb  ®R 

Diameter  of  the  duct,  ft 

Eccentricity 

(e/R) 

Length  of  the  duct,  ft 

Nondimension  length  of  the  duct,  (L/D) 

Mach  number 
(Uh/Ul) 

Mass  flow  rate,  Ib/sec 

* • 

Bypass  ratio,  {Ta.2/'^-^) 

Pressure 

Radius  of  the  duct,  ft 
Radial  coordinate 
Radius  of  the  free  jet 

Radius  of  the  primary  ejector  nozzle,  ft 

(y/R) 

Temperature,  "R 
Total  temperature,  *R 
(T  - Tj) 

Axial  velocity  component,  ft/sec 


iii 


AU 


(U  - U^) 

Axial  coordinate 
Transverse  coordinate 
Area  ratio,  (Aj^/A2) 

(Tm/T2) 

Temperature  ratio,  (T2/Tj^) 
Density,  Ib/ft^ 

(y/r) 

(R/r) 

3.14159 

Ratio  of  specific  heats 


Subscrii 


1 Primary  stream 

2 Secondary  stream 

3 Completely  uniform  mixed  flow 

h Nominal  wake  conditions 


m At  the > axis 

s Static 

w Wall  of  the  duct 


1 . INTRODUCTION 


Previous  analysis  has  shown  that  at  medium  and  low  flight 
Mach  numbers  an  exhaust  gas  jet  of  relatively  small  mass  and 
high  velocity  is  an  inefficient  method  of  producing  thrust. 

This  is  due  to  the  high  energy  losses  at  the  nozzle  exit  and 
the  poor  use  of  the  kinetic  energy  of  the  exhaust  jet  for  the 
production  of  thrust.  Two  other  important  considerations  for 
current  engines  are  good  fuel  economy  and  a low  noise  level. 

In  order  to  accomplish  these  objectives  of  maximizing  thrust, 
lowering  fuel  consumption  and  decreasing  noise,  a mixed  flow 
turbofan  engine  is  the  optimal  candidate. 

A high  airflow  is  necessary  to  increase  thrust  while  a 
high  jet  velocity  will  increase  thrust  but  decrease  propulsive 
efficiency.  The  turbofan  engine  produces  a compromise  between 
maximum  airflow  and  maximum  jet  exhaust  velocity.  Theoretical 
studies  have  shown  that  the  thrust  of  a turbofan  engine  can  be 
increased  by  mixing  the  hot  primary  jet  with  the  colder 
secondary  airflow.  Dr.  Hartmann  [1]  confirmed  these  theories 
with  his  experimental  investigations  at  a bypass  ratio  of  3.2. 

He  reported  a thrust  gain,  with  equal  total  pressure  in  the 
primary  and  secondary  streams  before  mixing,  for  mixing 
chamber  inlet  Mach  numbers  of  0.37  and  below.  The  amount  of 
thrust  gain  was  a function  of  the  stream  temperature  and  the 
degree  of  mixing. 

Therefore  it  is  important  to  optimize  the  amount  of  mixing 
in  a turbofan  nozzle.  In  order  to  optimize  this  mixing,  it  is 
necessary  to  better  understand  the  mechanics  of  turbulent  jet 
mixing.  Although  an  actual  turbofan  engine  may  have  turbine 
exit  swirl  and  an  augmentor  flame  holder  in  the  flow  path,  a 
simplified  analysis,  which  is  both  reliable  and  accurate,  is 
necessary  before  the  more  complex  problem  can  be  solved.  To 
satisfy  this  need  for  a tenable  turbulent  mixing  theory,  several 
people  have  made  analytical  studies  of  mixing  in  recent  years. 
Several  of  these  studies  will  be  discussed  in  this  paper.  The 
theory  proposed  by  Khanna  and  Tabakoff  [2] , however,  seemed  to 
provide  a fairly  simple  yet  promising  approach  to  the  turbulent 
mixing  problem.  The  theory  lacked  an  experimental  investigation 
of  the  predicted  results.  Therefore,  the  current  study  was 
performed  to  experimentally  verify  the  accuracy  of  the  theoretical 
analysis  of  turbulent  jet  mixing  between  two  concentric  air 
streams  in  a constant  area  duct.  In  particular  the  non-isoenergetic 
(different  stream  temperatures)  case  was  ex2unined.  An  eccentric 
mixing  configuration  was  used  as  an  additional  experimental 
study.  Non-isoenergetic  turbulent  jet  mixing  experiments  were 
conducted,  where  the  measured  temperatures  were  compared  with  a 
slightly  modified  version  of  Khanna  and  Tabakoff 's  theory.  These 
two  non-isoenergetic  mixing  studies  were  believed  to  have 
resulted  in  a better  understanding  of  turbulent  mixing  in  a 
turbofan  engine. 
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2.  LITERATURE  REVIEW  OF  RELATED  THEORIES  AND  EXPERIMENTS 


The  main  item  which  becomes  apparent  when  doing  a literature 
survey  is  the  lack  of  information  concerning  turbulent  ducted 
mixing  for  the  non-isoenergetic  case.  The  case  of  free  turbulent 
mixing  has  been  analyzed  by  many  authors,  such  as  Refs.  [3] , 

[4],  [5],  [6],  [7],  [8]  and  (9],  and  will  not  be  reiterated 

here.  The  case  of  ducted  isoenergetic  turbulent  mixing  is 
examined  in  several  works,  such  as  Refs.  [10],  [11],  [12],  [13], 
[14],  [15]  and  [16].  The  best  to  date  appears  to  be  that  of 
Razinsky  [17] . The  objective  of  his  investigation  was  to 
determine  the  velocity  field,  the  shear  stress,  and  the  pressure 
distribution  from  the  well  defined  step  change  in  the  velocity 
profile  at  the  mixing  tube  entrance  to  the  limiting  condition  of 
fully  developed  flow.  He  integrates  the  equations  of  motion  in 
the  radial  direction  for  preselected  forms  of  the  velocity 
profile.  The  assumptions  used  in  this  method  are  the  mean  radial 
velocity  is  small  compared  to  the  mean  axial  velocity,  gradients 
in  the  axial  direction  are  small  relative  to  the  radial  gradients, 
and  the  viscous  stress  is  small  compared  to  the  turbulent  shear 
stress  except  near  the  wall.  He  performed  an  experimental 
investigation  of  his  theory  which  generally  produced  good  agree- 
ment between  the  theory  and  the  experimental  data.  The  inclusion 
of  the  wall  boundary  layer  in  the  theoretical  model  yielded  a 
significantly  improved  prediction  of  the  pressure  distribution 
compared  to  previous  theories . 

Of  course  the  main  reference  for  this  paper.  Ref.  [2] , 
considered  both  isoenergetic  and  non-isoenergetic  turbulent 
mixing  from  an  analytical  viewpoint.  An  experimental  investigation 
of  the  isoenergetic  case  was  conducted  in  Ref.  [18]  . This  study 
agreed  with  the  theory  of  Ref.  [2]  for  the  case  where  the 
primary  stream  velocity  was  greater  than  the  secondary  stream 
velocity.  In  order  to  obtain  this  agreement,  however,  the 
coefficient  of  turbulence  (c)  had  to  be  modified.  Reference  [2] 
used  a value  of  c equal  to  0.2.  This  value  was  recommended  in 
Ref.  [3].  It  was  experimentally  determined  for  an  axially 
symmetric  free  jet.  For  the  case  of  confined  jets.  Ref.  [18] 
experimentally  determined  a value  of  c equal  to  0.7.  The 
confining  walls  require  the  mixing  to  have  this  new  constant  of 
turbulence,  which  is  completely  different  than  that  for  the 
free  jet  case.  In  order  to  provide  additional  understanding  of 
turbulent  mixing  an  experimental  Investigation  of  eccentric 
isoenergetic  ducted  mixing  was  conducted  in  Ref.  [19]  and  was 
compared  with  a slightly  modified  version  of  the  theory  of 
Ref.  [2].  This  work  was  expanded  by  the  non-isoenergetic  mixing 
experiments  discussed  in  this  paper. 
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3. 


THEORETICAL  ANALYSIS 


Discussion 


A review  of  the  theoretical  analysis  described  by  Khanna 
and  Tabakoff  [2]  is  appropriate  since  it  provided  the  basis  upon 
which  this  study  expanded.  He  found  that,  for  free  turbulent 
jets,  the  differential  equations  were  favorable  to  an  analytical 
solution  and  the  experimental  behavior  was  quite  well  known. 

Much  less  well  known,  however,  was  the  behavior  of  a jet  immersed 
in  a secondary  stream  with  confining  walls  present.  Little 
experimental  information  on  turbulent  shear  stress  was  available 
for  jet  flows  other  than  free  jets.  Extensive  measurements  on 
free  jets  had  shown  that  this  type  of  flow  was  closely  self- 
preserving, up  to  a certain  distance  downstream  of  the  nozzle 
exit.  Using  the  ideas  of  Abramovich  [3]  that  hypothesized  that 
the  velocity  or  temperature  profile  in  a duct  would  be  the  portion 
of  a free  jet  profile  located  between  the  duct  walls,  he  developed 
the  following  analysis  for  turbulent  jet  mixing  between  two 
subsonic  streams  in  an  axisymmetric  constant  area  duct.  His 
analysis  included  both  the  constant  temperature  (isoenergetic) 
mixing  case  and  the  different  temperature  for  each  stream  (non- 
isoenergetic)  mixing  case. 

Figure  1 illustrates  the  schematics  of  the  ducted  mixing 
system.  It  is  divided  into  three  flow  regimes  which  are  defined 
as  follows: 

1.  First  regime  or  initial  region.  In  this  region  turbulent 
mixing  occurs  between  the  secondary  flow  and  the  core  of  inviscid 
primary  flow.  When  the  inner  boundary  of  this  mixing  region 
intersects  with  the  axis  of  the  jet,  the  region  is  ended. 

2.  Second  regime  or  transitional  region.  In  this  region 
the  inviscid  core  of  primary  flow  has  been  dissipated,  but  a 
region  of  inviscid  secondary  flow  exists  near  the  mixing  chamber 
wall. 


3.  Third  regime  or  main  region.  In  this  region  the  mixing 
layer  has  spread  to  the  wall  and  no  portion  of  undisturbed 
secondary  flow  exists. 

It  is  generally  considered  that  the  first  and  second  regimes 
make  up  the  initial  region  and  that  the  third  regime  is  consi- 
dered the  main  region.  The  theory  used  here  applies  only  to  the 
main  region.  Other  assumptions  used  by  this  theory  are:  viscous 

effects  at  the  wall  are  negligible;  the  static  pressure  at  any 
cross-section  is  constant;  mixing  obeys  the  perfect  gas  laws; 
and  the  duct  walls  are  adiabatic. 
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Derivation  of  Equations 

Since  the  general  analytical  problem  is  complex,  the  theory 
first  addressed  the  isoenergetic  case.  Later  the  appropriate 
changes  and  additions  were  made  for  the  non-isoenergetic  case. 
The  fundamental  equations  governing  the  flow  in  ducts  are; 

1.  Conservation  of  Mass.  For  one-dimensional  steady  flow 
the  continuity  equation  can  be  written  as 

3(pU)  ^ 0 

3x 

or  pAU  = constant  = m 
applying  this  equation  to  our  model,  we  have 

Pi  Ai  Ui  + P2  A2  U2  = P3  A3  U3  (1) 


2.  Conservation  of  Momentum:  Neglecting  the  viscous 

forces  acting  on  the  walls  of  the  mixing  chamber,  we  can  write 
the  momentum  balance  equation  as 

(P3  - PjjAj  - P3  Vl  A3  + pj  Aj  - pj  A3  (2) 


3.  Conservation  of  Energy:  Assuming  the  specific  heat  at 

constant  pressure  to  be  constant,  independent  of  the  temperature, 
the  equation  of  conservation  of  energy  is 


m. 


P 


'^p  ^t. 


= m. 


(3) 


4.  Equation  of  State;  For  a perfect  gas,  the  equation  of 
state  can  be  written  as 

p = pRT  (4) 

Equations  (1)  to  (4)  determine  the  flow  parameters  after  complete 
mixing,  which  occurs  at  an  infinite  distance  from  the  initial 
cross-section  of  the  mixing  chamber.  The  calculations  of  flow 
parameters  for  a short  mixing  chamber,  in  which  we  have  to  take  into 
account  the  incomplete  mixing,  and  the  determination  of  the 
optimum  length  of  the  mixing  chamber,  requires  the  knowledge  of 
the  laws  of  mixing  for  coflowing  streams  along  the  length  of  the 
mixing  chamber. 


4 


In  order  to  determine  the  variation  of  stream  velocity 
and  temperature  along  the  radius  of  the  mixing  chamber  at  any 
cross-section,  we  shall  make  use  of  the  hypothesis  of  universal 
ejection  characteristics  developed  by  Abramovich  [3] . This 
hypothesis  states  that  the  nondimensional  velocity  and  temperature 
profiles  at  any  cross-section  of  a turbulent  jet  are  universal 
functions  of  a parameter  5 = (y/r) , regardless  of  the  external 
conditions  of  its  development.  In  the  C term,  y is  the 
distance  from  the  axis  in  the  transverse  direction  and  r is  the 
radius  of  the  free  jet.  This  hypothesis  allows  us  to  use  the 
results  of  free-jet  mixing  theories  to  determine  the  velocity 
and  temperature  profiles  for  turbulent  mixing  of  two  streams 
in  a duct.  Mathematically,  we  can  write  the  universal  functions 
for  nondimensional  excess  velocity  and  excess  temperature 
profiles  as; 


AU 

AUm 


(1  - 


(5) 


AT 

ATm 


(1  - 5^*^) 


(6) 


For  axisymmetric  turbulent  mixing  of  free  jets,  C = (y/r).  But 
for  mixing  of  two  streams  in  a constant  area  mixing  chamber, 
the  maximum  value  of  y is  equal  to  R,  the  radius  of  the  mixing 
chamber.  The  corresponding  value  of  C is  denoted  by  = (R/r)  , 

which  is  always  less  than  one  in  the  main  region  of  the  mixing 
chamber.  Figure  2 better  defines  the  terms  used  in  these 
equations.  This  figure  illustrates  the  basic  concept  used  to 
find  the  velocity  and  temperature  profile  inside  the  duct. 

The  excess  velocity  at  the  axis  AUm,  and  the  excess  velocity  AU, 
are  defined  by 


AUm  = Urn  - Uh 
AU  = U - Uh 


(7) 


where  Urn  = axial  velocity  at  the  axis 
Uh  = nominal  wake  velocity 
U = axial  velocity  at  any  value  of  y 

The  nominal  wake  velocity,  Uh,  can  be  obtained  from  Bernoulli's 
equation  if  it  is  assumed  that  the  main  change  in  the  pressure 
occurs  in  the  initial  region,  and  that  the  longitudinal  pressure 
gradient  in  the  main  region  is  small,  i.e.,  a constant  pressure 
is  assumed  for  the  whole  main  region.  Therefore, 


Uh^  = ^2  " f ^^3  " ^2^ 

where  p,  is  the  pressure  after  completely  uniform  mixed  flow  has 
been  attained. 


I 
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Non-Isoenergetic  Mixing 

The  analysis  of  the  previous  section  will  now  be  extended 
to  include  the  case  of  an  initial  temperature  difference  between 
the  two  mixing  airflows.  This  was  the  primary  case  of  concern 
to  this  author,  since  the  purpose  of  the  experiments  was  to 
measure  temperatures,  to  verify  the  theoretical  analysis  and  to 
predict  the  behavior  of  non-isoenergetic  turbulent  mixing  in 
a duct. 

One  can  define  a criterion  for  the  degree  of  mixing  as  the 
ratio  of  the  total  energy  actually  absorbed  by  the  lower  energy 
stream  to  the  maximum  absorbable  total  energy,  which  corresponds 
to  perfect  mixing  of  the  streams  with  complete  equalization  of 
states  of  the  mixing  gases.  This  essentially  is  a measure  of  the 
degree  of  uniformity  of  the  temperature  profile  in  the  mixing 
duct.  The  variation  of  temperature  along  the  radius  of  the  mixing 
duct  at  any  cross-section  will  be  determined  by  the  universal 
function  of  the  parameter  (C)  as  stated  before.  To  determine 
the  variation  of  temperature  along  the  length  of  the  duct,  we 
shall  write  the  integral  relation  which  expresses  the  law  of 
conservation  of  mass  as 


I 


I 


j 


/ pU  dA  = Aj 

o 


Now  adding  and  subtracting  the  quantity 
^3 

/ pUh  dA 
o 

from  Equation  (9) , we  have 


Since 


"3  *^3 

f p (U  - Uh)dA  + / pUh  dA 
o o 


^3  ^3  ^3 


= irR 


and 


dA  = 2ity  dy 


(9) 


(10) 
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Therefore,  Equation  (10)  becomes 


i 


R R 2 

2 I pAUy  dy  + 2Uh  /pydy  = p^U-R 

o o 

For  approximately  constant  pressure  mixing,  using  the 
equation  of  state.  Equation  (4),  we  have 


(11) 


= E_ 

RT 


(12) 


Substituting  Equation  (4)  into  (11) , we  have 


Or 


Or 


2 / (AU/T)y  dy  + 2Uh  / (y/T)dy  = R /T- 

o o 


(13) 


2 iUm  / + 20h  / = U3  ^(Tm/Tj) 

(14) 


• • 


Now  (pj^/p2)  = 0 = " “ (m2/ni2^)  and  denoting  1))  = (Tm/T2) 

we  have 


Tm  _ Tm  "^2  ^1  _ ^ „ ,n  + 1 ^ 

T3  ~ T2  Tj^  Tj  “ ’’  ® 'ne  + 


(15) 


Using  Equation  C6)  and  the  expression  for  (fi  we  have 

CT/Tml  = II  - C-^)  5^*^] 

4 


(16) 


Substituting  Equation  (16) , (15)  and  (5)  into  Equation 
(14)  we  have 


,1.5»2 


, (l-C*’-")‘  6dC  . Uh  , 

^ ITT"  + — ^ 


AOm 


CdC  „ 1 ytt+1  \ A / j , 1 » 

T3  I 


(17) 


i 


Integration  of  Equation  (17)  produces  the  following  equation, 

3 

4 ,o(l+ne)-mCa+l)  1 r k B -2  ,_3  . B „ 

^ ^0~(n+T5  ^ ‘T-  + —4"  F ^1 


, b/7  Tin  , 4m(a+l)  rv  . S „ , b/T  tt, 

T8— aeTn+ry  ^^k  F ^1  -T8~’ 


1-B' 


(18) 


where : 


(*/^_l)  V3 

(X.  - B)^ 

In  — 5 *—  - 2/3  tan 

(Xj^  + BXj^  + B ) 

= 1.0  - 1.143  + 0.4 


2X.  + B 

(_JS 

Bv'I 


) 


a = (Aj^/A2) 

m = (Uh/Uj^)  (19) 


Equation  (18)  implicitly  determines  the  quantity  (Tm/T2) , 
where  Tm  is  the  temperature  on  the  axis  of  the  mixing  duct 
at  any  cross-section.  As  can  be  seen  from  Equation  (18)  , the 
quantity  ^ appears  on  both  sides  of  the  equation  and  therefore 
it  is  not  possible  to  determine  the  temperature  on  the  axis  of 
the  mixing  duct  for  any  specified  location  of  the  cross-section. 
Instead,  a trial  and  error  scheme  is  required  to  determine  the 
temperature  distribution  along  the  length  of  the  mixing  duct. 

The  details  of  the  iterative  scheme  employed  to  solve  Equation 
(18)  on  the  CDC  6600  computer  are  given  in  Appendix  A. 


4.  TEST  CONFIGURATION  AKD  DESCRIPTION 
Apparatus  and  Test  Facilities 

A.  Description  of  the  test  stand: 

A schematic  diagraun  of  the  mixing  tunnel  is  shown  in  Figure  3. 
Air  is  supplied  from  an  air  pressure  tank  at  200  psia.  The  air 
passes  through  a pressure  regulator  in  order  to  reduce  the  pressure 
of  the  flow  to  the  desired  value.  As  a rule,  the  pressure  after 
the  pressure  regulator  was  held  constant  at  100  psia  during  each 
run,  and  was  checked  during  the  runs  with  the  flow  control  panel. 
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The  total  air  flow  is  measured  with  an  orifice  meter.  The  main 
flow  is  then  divided  into  two  flows,  the  primary  and  the  secondary 
flows.  The  amount  of  mass  flow  in  each  is  controlled  with  valves. 

The  primary  flow  is  routed  through  a pipe  to  the  combustion 
chamber,  which  was  used  to  heat  the  primary  airflow.  With  the 
use  of  several  different  fuel  flow  nozzles  and  a fuel  flow  control 
panel,  it  was  possible  to  regulate  the  primary  flow  temperature. 

The  heated  flow  then  passed  through  an  insulated  pipe  until  it 
arrived  at  the  entrance  to  the  mixing  tunnel.  The  secondary  air 
is  measured  with  a second  orifice  meter,  then  it  enters  a 
settling  chamber,  in  order  to  change  its  direction. 

The  primary  flow  passes  through  the  settling  chamber  inside 
a two-inch  diameter  pipe  which  could  easily  be  exchanged  with 
another  of  diameter  1.4  inches.  The  primary  flow  pipe  has  a 
sharp  edge  so  that  the  thickness  at  the  initial  mixing  cross- 
section  is  practically  zero.  The  secondary  flow  is  then  allowed 
to  flow  co-axially  with  the  primary  flow  inside  the  mixing  duct 
which  extends  five  feet  downstream.  A movable  rake  was  provided 
in  the  mixing  duct.  The  mixing  duct  inside  diameter  was  four 
inches. 

The  flow  control  panel  is  connected  to  pressure  lines  upstream 
of  two  orifice  meters  and  is  used  to  assure  that  the  air  pressure 
after  the  regulator  is  kept  constant.  A switch  is  connected  to 
the  flow  control  panel  which  could  open  the  air  line  instantaneously. 


B.  Instrumentation 

The  air  used  during  the  test  was  supplied  from  a storage 
tank.  During  each  run,  the  pressure  inside  the  tank  dropped 
continuously.  If  the  pressure  inside  the  tank  began  to  drop 
below  the  pressure  set  up  by  the  regulator,  the  mass  flow  of 
air  began  to  change,  and  the  conditions  specified  for  a run 
changed.  Therefore,  the  duration  of  the  test  was  limited.  In 
addition,  it  required  between  ten  and  twenty  minutes  to  regulate 
and  then  stabilize  the  desired  temperature  of  the  primary  flow. 
Monitoring  of  this  temperature  was  also  required  during  the  run 
in  order  to  maintain  a constant  temperature.  Thus,  it  was 
essential  to  have  a system  of  pressure  and  temperature  probes 
that  could  quickly  measure  the  radial  distribution  at  each 
cross-section.  Therefore,  a temperature/pressure  rake  was  built 
for  the  measurements . 

A schematic  drawing  of  the  rake  used  is  shown  in  Figure  4 . 
Seven  stainless  steel  tubes,  each  with  an  outside  diameter  of 
0.13  inches,  were  connected  to  the  rake.  Five  of  the  tubes 
were  closed  at  the  end  with  a copper-constantin  thermocouple. 

The  other  two  tubes  were  used  to  measure  the  total  pressure. 

The  tubes  were  equally  spaced  from  each  other  at  multiples  of 
0.5  Inches  from  the  duct  centerline.  Each  pipe  was  secured  at 
its  proper  distance  inside  the  wedge-shaped  rake.  The  rake 
shaft  was  installed  to  insure  that  it  moved  exactly  along  the 
mixing  duct  axis.  The  rake  was  connected  to  a shaft  inside  of 
which  passed  the  pressure  probes  and  thermocouple  wires.  They 
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ran  from  the  rake  to  the  appropriate  data  recorder  instru- 
ment. The  pressure  probes  were  connected  to  a multi-tube 
manometer  which  also  recorded  the  static  pressure.  The 
thermocouple  wires  were  attached  to  a temperature  recorder. 

This  instrument  converted  the  millivolt  signals  from  the  wires 
into  the  corresponding  temperatures.  These  temperatures  were 
then  plotted  on  grid  paper  in  a different  color  and  number  for 
each  thermocouple.  Thus,  the  temperatures  could  be  read 
directly  from  the  recorder  in  degrees  Fahrenheit.  The  entire 
rake  system  slid  along  a scale  so  that  an  accurate  position  of 
the  rake  inside  the  mixing  duct  was  known. 

This  test  apparatus  and  instrumentation  for  the  eccentric 
case  was  identical  to  the  concentric  case  with  the  exception 
of  the  position  of  the  primary  flow  tube.  Both  the  2.0  inch 
and  the  1.4  inch  diameter  primary  flow  tubes  were  positioned 
inside  of  the  4.0  inch  diameter  duct  in  such  a manner  that  their 
centerlines  were  displaced  one-half  inch  below  the  duct  center- 
line.  In  the  concentric  case  the  centerline  was  identical  for 
both  tubes.  Figure  5 illustrates  the  eccentric  configuration. 

Testing  Procedure  and  Accuracy 

During  each  run  the  rake  was  set  at  a certain  position  and 
the  multitube  manometer  was  manually  read.  The  temperature 
recorder  worked  automatically  for  the  complete  run.  Readings 
were  taken  at  only  five  axial  locations  due  to  the  run  time 
limitation.  The  five  locations  corresponded  to  five  cross- 
sections  located  downstream  of  the  initial  mixing  cross-section 
at  distances  of  0.2,  3,  5,  8,  and  12  times  the  internal  diameter 
of  the  mixing  duct.  This  procedure  permitted  fairly  quick 
measurements  with  a temperature  error  of  approximately  plus 
or  minus  five  degrees  Fahrenheit  and  a pressure  error  not 
greater  than  five  percent  for  the  case  where  the  velocity  was 
between  200  and  400  ft/sec. 

The  velocity  was  readily  calculated  from  the  formula 

2 2 Pt  ^ 

M = ^ [(^)  - 1]  (20) 

Since  both  the  total  and  static  pressure  and  temperature  at  a 
given  point  were  known. 

Experimental  Results 

Three  different  test  configurations  were  employed  during  the 
total  experiment.  These  were:  (1)  Area  Ratio  = 3.0  and  U1  > U2; 

(2)  Area  Ratio  =*  3.0  and  Ul  < U2;  (3)  Area  Ratio  = 7.16  and 

U1  > U2. 


10 


(A)  Concentric  Configuration  - The  experimental  temperature 
profiles  are  plotted  in  Figure  6 for  the  Area  Ratio  = 3.0, 

U1/U2  = 2.4  and  T1/T2  = 2.0  case  at  the  five  different  cross- 
sections.  Results  for  other  temperature  ratios  can  be  found 
in  Ref.  (24].  A fully  developed  temperature  profile  was  obtained 
for  most  of  the  results  at  an  L/D  equal  to  twelve  (additional 
experimental  results  are  given  in  Ref.  [24]).  This  means  that 
the  nonuniformity  of  the  temperature  profile  has  disappeared 
due  to  the  high  degree  of  mixing  in  the  duct.  From  the  data 
obtained  we  find  that  the  larger  the  initial  temperature  ratio, 
the  longer  it  takes  for  the  flow  to  become  completely  mixed. 

It  was  also  observed  that  the  temperature  nonuniformity  disappears 
more  rapidly  than  the  velocity  nonuniformity.  Velocity  profiles 
were  not  plotted  since  they  are  given  in  Ref.  (24]. 

Figure  7 is  a plot  of  the  experimental  temperature  profiles 
for  the  U1  > U2  case  where  the  Area  Ratio  = 7.16.  This  area 
ratio  was  obtained  by  using  a 1.4  inch  dieuneter  pipe  inside 
the  4.0  inch  diameter  duct.  Again  the  data  taken  at  the  five 
axial  cross-sections  is  plotted.  A fully  developed  temperature 
profile  is  again  obtained  at  a distance  equal  to  twelve  times 
the  diameter  of  the  mixing  duct  (L/D  = 12.0).  The  rate  of  decay 
of  the  temperature  of  the  primary  flow  is  higher  than  for  the 
previous  case  where  the  area  ratio  was  equal  to  three.  In  both 
the  U1  greater  than  U2  cases,  it  was  expected  that  the  greater 
the  velocity  ratio  or  the  temperature  ratio,  the  greater  the 
degree  of  turbulence,  which  would  decrease  the  profile  development 
length. 

Figure  8 is  a plot  of  the  experimental  temperature  profiles 
for  the  U1  less  than  U2  case  where  the  area  ratio  equals  3.0, 
the  velocity  ratio  is  0.59  and  the  temperature  ratio  is  1.65. 

The  data  shows  that  the  final  temperature  is  heavily  influenced 
by  the  secondary  stream  temperature.  This  is  due  to  the  fact 
the  bypass  ratio  is  high  and  the  primary  flow  is  rapidly  damped 
out.  The  low  temperature  ratios  indicated  more  rapid  mixing  in 
the  U1  greater  than  U2  cases,  and  this  is  also  true  for  the 
U1  less  than  U2  case. 

(B)  Eccentric  Configuration  - The  first  case  is  again  the 
U1  greater  than  U2  configuration  at  an  area  ratio  equal  to  3.0. 
Figure  9 represents  the  experimental  temperature  profiles  at 
the  five  axial  locations  along  the  duct  which  were  measured  for 
this  case.  Figure  10  shows  a configuration  where  U1  is  greater 
than  U2,  but  the  area  ratio  is  equal  to  7.16.  The  temperature 
and  velocity  ratio  values  are  a repeat  of  the  same  values  used 
for  concentric  mixing.  Figure  11  represents  the  U1  less  than 
U2  configuration  at  an  area  ratio  of  3.0.  The  same  velocity 
ratios  as  the  concentric  experiments  were  used. 

From  the  temperature  curves  presented  in  Figures  9,  10  and 
11  it  is  apparent  that  the  temperature  profiles  are  always  un- 
symmetrical  about  the  axis  of  the  mixing  duct.  At  the  cross- 
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sections  near  the  initial  cross-section,  the  maximum  temperature 
is  at  the  axis  of  the  primary  flow.  But  as  one  proceeds  down- 
stream the  peak  temperature  tends  to  move  toward  the  wall  on  the 
side  of  the  eccentricity.  Additionally,  the  temperature  non- 
uniformity does  not  decay  with  the  same  rate  as  in  the  case  of 
concentric  mixing.  The  temperature  decays  more  rapidly  with 
increasing  area  ratio  and  when  U1  is  less  than  U2,  this  agrees 
with  the  expected  result.  All  three  configurations  behave  in 
the  same  general  fashion  with  regard  to  temperature  profile 
decay.  The  profile  remains  unsymmetrical  even  at  the  final 
cross-section  measured  (i.e.,  L/D  = 12.0).  It  is  seen  that 
eccentricity  increases  the  axial  distance  required  for  a fully 
developed  profile  to  be  obtained.  The  rate  of  decay  of  the 
temperature  profile  does  not  appear  to  have  a direct  relation- 
ship to  either  the  prim.ary  temperature  magnitude  or  the 
temperature  ratio.  The  rate  of  decay  appears  to  be  fairly 
consistent  with  respect  to  these  variables.  This  is  again 
different  than  the  velocity  mixing  studies,  where  the  rate 
of  decay  of  the  velocity  at  the  primary  flow  axis  increased 
with  both  increasing  velocity  ratio  and  velocity  magnitude. 


5.  RESULTS  AND  DISCUSSION 

A preliminary  examination  of  the  test  data  indicated  that 
the  temperature  mixing  was  much  more  rapid  than  the  theory 
predicted.  Therefore,  based  on  the  experience  of  Ref.  [18] 
where  a new  value  of  the  constant  of  turbulence  (c)  for  the 
main  region  was  experimentally  obtained,  a computer  check  was 
made  on  other  values  of  c.  It  was  found  that  increasing  the 
value  of  c,  even  to  a value  of  20,  did  not  significantly  improve 
the  agreement  between  the  theoretical  and  experimental  data. 
Another  approach  was  deemed  necessary.  Following  the  background 
and  reasoning  of  Ref.  [18]  , the  values  of  the  temperature  at 
the  axis  (Tm)  and  the  temperature  near  the  wall  outside  the 
boundary  layer  (Tw)  were  plotted  for  various  values  of  L/D 
along  the  mixing  duct.  Instead  of  changing  the  value  of  c, 
however,  the  value  of  Cj^,  the  constant  of  turbulence  in  the 

initial  region,  was  varied  and  the  resulting  theoretical  values 
of  Tm  and  Tw  were  plotted.  Since  it  was  apparent  that  the 
initial  temperature  difference  of  the  two  streams  caused  the 
mixing  to  be  accelerated,  it  was  reasonable  to  assume  that  the 
accelerated  mixing  would  also  affect  the  mixing  in  the  initial 
region  and  hence  shorten  the  axial  length  of  the  initial  region. 
Various  values  of  Cjj  were  plotted  against  the  experimental  data 

in  Figure  12.  A value  of  Cj^  =-0.36  presented  the  best  agreement 

between  the  theory  and  the  experimental  data.  The  case  where 
U1  was  less  than  U2  was  not  considered  significant  due  to  the 
negative  results  of  [18]  and  an  exaunination  of  our  experimental 
test  results,  which  showed  that  good  agreement  was  not  obtained. 

With  the  value  of  c^  » 0.36  and  the  value  of  c “ 0.7  the 
computer  program  was  run  to  obtain  the  theoretical  curves  used 


to  compare  with  the  experimental  data. 


After  having  examined  the  experimental  temperature  profile 
plots  and  the  theoretical  versus  experimental  temperature  plots, 

. several  observations  can  be  made.  The  accuracy  of  the  temperature 

1 measurement  system  should  be  kept  in  mind  when  comparing  various 

temperature  values.  Within  these  restraints  the  data  agrees 
reasonably  well  with  the  experiments  for  the  cases  where  U1  is 
greater  than  U2.  There  can  be  seen  a definite  trend  for  the 
experimental  data  to  have  a lower  temperature  than  the  theoretical 
data  at  higher  L/D  values.  The  cause  of  this  observed 
phenomenon  is  believed  to  be  the  fact  that  heat  loss  occurred 
through  the  duct  walls  during  the  tests.  Thus  the  adiabatic 
wall  assumption  is  probably  not  valid  for  these  experiments. 

This  phenomenon  will  be  examined  further  in  the  subsequent 
sections . 

Theoretical  versus  Experimental  mixing  comparison  for  the 
concentric  case  - As  was  previously  discussed,  a new  value  of 
the  constant  of  turbulence  for  the  initial  region  was  required 
to  account  for  the  increased  degree  of  mixing  which  was  exhibited 

by  the  experimental  data.  This  new  value  helped  in  the  agree- 

ment of  the  theoretical  and  experimental  data  for  the  cases 
where  U1  was  greater  than  U2.  it  did  not  help  in  the  U1  less 
than  U2  cases.  Using  the  new  value  of  the  constant  of 
turbulence  (Ch  = 0.36)  a comparison  between  theory  and  experi- 
ment was  made  with  plots  for  the  following  cases.  It  should 
be  noted  that  T1  is  greater  than  T2  for  all  mixing  cases. 

Case  I:  U1  > U2,  Area  Ratio  = 3.0 

Figure  13  illustrates  this  case  for  velocity  ratio  of  2.4 
approximately.  Each  curve  on  the  figure  represents  a temperature 
profile  at  the  given  cross-section  for  different  values  of 
temperature  ratio.  The  rapid  degree  of  the  disappearance  of 
the  initial  temperature  nonuniformity  is  readily  apparent. 

At  large  values  of  L/D  (Ref.  [24]),  the  experimental  temperature 
profile  has  a lower  mean  temperature  than  the  theoretical 
temperature  profile.  This  phenomenon  was  due  to  the  heat  loss 
through  the  wall  of  the  duct.  In  general,  the  agreement  between 
theory  and  experiment  is  fairly  good.  The  velocity  ratio  does 
not  seem  to  have  much  effect  on  the  degree  of  mixing  within  the 
accuracy  limits  of  the  data.  As  was  previously  mentioned  the 
temperature  profile  becomes  fully  developed  more  quickly  for  a 
lower  temperature  ratio. 

Case  II:  U1  > U2,  Area  Ratio  = 7.16 

I 

■I  This  case  is  illustrated  by  Figure  14.  Again  a plot  is 

made  for  a given  velocity  ratio  and  a given  L/D  with  curves  of 
two  different  temperature  ratios.  Because  of  the  higher  area 
ratio,  the  initial  region  length  is  decreased.  Case  II  figure 
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shows  a lesser  effect  of  heat  loss.  Only  a small  decrease  in 
the  experimental  temperature  data  is  seen  at  an  L/D  equal  to 
twelve  (Ref.  [24]).  This  is  the  expected  result  since  the  mass 
flow  of  cold  secondary  air  compared  to  the  mass  flow  of  hot 
primary  air  is  significantly  larger  than  for  Case  I.  The  lower 
average  temperature  in  the  partially  mixed  flow  contributes  to 
much  less  heat  loss  through  the  duct  wall.  The  agreement 
between  the  Case  II  experimental  and  theoretical  temperature 
profiles  is  quite  good. 


6.  CONCLUSIONS 

An  experimental  investigation  of  non-isoenergetic  turbulent 
jet  mixing  between  two  streams  in  an  axisymmetric  duct  has  been 
conducted.  Experimental  temperature  profiles  have  been  measured 
at  several  axial  locations  in  the  duct.  They  have  been  compared 
with  the  theoretical  temperature  profiles  determined  by  the 
analysis  described  in  this  paper.  Several  conclusions  can 
be  made  from  these  comparisons. 

It  was  determined  that  the  theory  and  experimental  results 
agree  fairly  well  for  the  concentric  case  where  the  primary 
stream  velocity  is  greater  than  the  secondary  streeim  velocity. 
This  is  true  for  both  area  ratios  tested.  The  velocity  ratio 
had  no  significant  effect  on  the  rate  of  temperature  mixing  of 
the  flow.  Also,  contrary  to  the  theory,  the  mixing  process 
was  not  accelerated  by  a larger  initial  temperature  difference 
between  the  flows.  The  opposite  effect  appeared  to  be  the  actual 
effect.  It  was  experimentally  verified,  however,  that  energy 
diffuses  more  rapidly  than  momentum.  The  rate  of  mixing 
determined  by  the  concentric  mixing  experiments  indicated  that 
the  mixing  was  more  rapid  than  the  theoretical  analysis  had 
predicted. 

The  theory  used  in  this  analysis  is  a fairly  simple  one 
which  produces  reasonably  good  results.  Eliminating  a couple 
of  the  simplifying  assumptions  would  provide  better  agreement 
with  the  experimental  data.  However,  it  is  not  thought  that  the 
increased  complexity  is  worth  the  additional  effort. 
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Inviscid  Secondary  Flow 


FIGURE  1 Schematic  of  Ducted  Mixing  System 


FIGURE  2 Schematic  Diagram  of  Confined  Jet  Mixing 


FIGURE  3 Mixing  Tunnel  Schematic 
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